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ABSTRACT 

This  paper  addresses  the  problem  to  define  a methodology  for 
the  analysis  of  the  performances  of  different  aircraft 
configurations  in  the  phase  of  conceptual  design.  The 
proposed  approach  is  based  on  a numerical  optimisation 
procedure  where  a scalar  objective  function,  the  take-off 
weight,  is  minimised.  The  optimisation  algorithm  has 
obviously  important  consequences  both  from  the  point  of  view 
of  the  computational  times  and  of  the  obtained  results . For 
this  reason  a preliminary  discussion  is  made  where  various 
different  methodologies  are  critically  compared.  Although  the 
best  compromise  between  different  approaches  is  probably 
given  by  an  integration  between  a genetic  algorithm 
approach  and  a classical  gradient  method,  in  this  phase  only 
the  latter  procedure  has  been  used  to  perform  the  simulations. 
The  methodology  takes  into  account  the  high  number  of 
geometrical  parameters  and  the  flight  mechanics 
requirements  involved  in  the  problem. 

A basic  example  is  described,  and  the  use  of  the  proposed 
methodology  to  investigate  the  effects  of  different  geometrical 
and  technological  parameters  is  discussed. 

NOMENCLATURE 

A Sectional  area  of  panel  structure,  m2 

AR  aspect  ratio 

b span,  m 

c chord  length,  m 

cD  global  drag  coefficient 

ci  sectional  lift  coefficient 

cL  global  lift  coefficient 

cma  mean  aerodynamic  chord  of  the  wing,  m 

E efficiency  (lift-to-drag  ratio) 

g acceleration  of  gravity,  m/s2 

H cruise  altitude,  m 

h-hn  distance  (non-dimentionalised  with  cma)  between 
aircraft  c.g.  and  neutral  point 
i stabiliser  angle,  deg 

L length,  m 

lt  distance  between  aerodynamic  centres  of  wing  and 
stabiliser,  m 


M Mach  number 

q dynamic  pressure,  N/m2 

R range  at  cruise  velocity,  km 

S surface  area,  m2 

Sjt  specific  fuel  consumption,  (N/hr)/N 

T vertical  gap,  m 

t/c  thickness  ratio 

V cruise  velocity,  m/s 

W weight,  N 

WS  wing  loading,  N/m2 

xcg  c.g.  position  of  the  residual  weight,  (%  Lflls) 
xs  position  of  stabiliser  root  leading  edge  from  fuselage 
nose,  m 

y spanwise  coordinate,  m 

a angle  of  attack,  deg 

F dihedral  angle,  deg 

<5  elevator  angle,  deg 

0 twist  angle,  positive  increases  the  tip  angle  of  attack, 
deg 

A taper  ratio 

A sweep  angle,  deg 

p air  density,  kg/m3 

pM  density  of  structural  material,  kg/m3 

oa  admissible  normal  stress  level,  N/m2 

0 diameter,  m 

Subscripts 

BL  bending  due  to  lift 

c cruise 

ea  elastic  axis 

/ fuel 

fus  fuselage 

M pitching  moment 

o take-off 

p payload 

r root 

ref  reference 

res  residual 

rib  wing  ribs 

5 stabiliser  surface 

SL  shear  due  to  lift 

t tip 

w wing 


Paper  presented  at  the  RTO  AVT  Symposium  on  “ Aerodynamic  Design  and  Optimisation  of  Flight  Vehicles  in  a 
Concurrent  Multi-Disciplinaiy  Environment ”,  held  in  Ottawa , Canada,  18-21  October  1999,  and  published  in  RTO  MP-35. 
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1.  INTRODUCTION 

During  the  phase  of  conceptual  design,  the  problem  of 
evaluating  the  performance  of  different  aircraft  configurations 
is  of  great  interest;  however,  difficulties  arise  due  to  the  high 
number  of  geometrical  parameters  involved,  which  are 
necessary  for  defining  such  a configuration.  A systematic 
analysis  taking  into  account  the  effects  of  all  these  parameters 
appears  to  be  difficult,  given  the  complexity  related  to  both 
aerodynamic  load  evaluation  and  the  assessment  of  flight 
mechanics  requirements.  Furthermore,  a sufficiently  precise 
estimation  of  the  weights  of  the  different  aircraft  components 
is  necessary.  However  this  estimation  is  strongly  related  to 
aircraft  lift  distribution,  hence  the  use  of  a satisfactory 
aerodynamic  code  is  mandatory. 

In  order  to  tackle  this  problem,  a direct  numerical 
optimisation  technique  may  be  employed,  giving  rise  to  an 
attractive  approach,  since  the  problem  may  be  addressed 
systematically,  and  the  designer  has  a great  flexibility  in  the 
choice  of  the  design  variables. 

In  the  analysis  through  direct  numerical  optimisation,  an 
aerodynamic  code  is  coupled  in  a loop  with  an  optimisation 
routine,  so  as  to  automatically  manage  the  values  of  the  design 
variables  - typically  concerning  geometry  modifications  - 
with  the  aim  of  minimising  a given  scalar  quantity  (objective 
function).  This  approach  is  extremely  flexible,  and  capable  of 
meeting  multidisciplinary  requirements.  The  choice  of  the 
design  variables  is  of  fundamental  importance,  as  they  define 
the  set  of  solutions  within  which  the  optimal  one  is  sought. 
Finally,  the  design  problem  formulation  requires  the  definition 
of  the  constraint  functions,  which  can  be  of  aerodynamic, 
flight  mechanics  or  geometrical  nature,  or  more  generally 
involve  any  quantity  that  may  be  computed  with  sufficient 
reliability. 

In  a previous  paper1  an  aerodynamic  optimisation 
procedure,  based  on  a potential  flow  model,  has  been 
developed  for  the  study  of  the  effects  of  the  geometric 
parameters  on  the  induced  drag.  The  obtained  code  proved  to 
be  quite  reliable  in  terms  of  efficiency  in  determining  the 
optimum  configuration  for  the  induced  drag,  but,  in  the 
version  presented  in  Ref,  1,  only  this  aerodynamic  aspect  was 
taken  into  account.  Consequently,  the  analysed  configurations 
did  not  meet  the  constraints  resulting  from  aircraft  flight 
mechanics,  nor  was  the  structural  weight  considered.  In  Ref.  2 
the  above  obstacles  were  overcame  by  including  the 
constraints  of  a realistic  configuration  in  the  code,  with  the 
assessment  of  trim  conditions,  static  stability  requirements  and 
with  the  addition  of  the  evaluation  of  airplane  structural 
weight,  based  on  a methodology  originally  developed  by 
Torenbeek.3  The  weight  was  included  with  a low  level  routine, 
the  aim  being  to  develop  a code  suitable  for  a conceptual 
design  approach.  This  simplification,  however,  should  not  be 
considered  a real  limitation  in  so  far  as  different  modules, 
which  make  up  the  code,  can  be  substituted  by  other  ones, 
with  more  sophisticated  mathematical  models,  without  the 
need  to  modify  the  problem  approach. 

This  paper  is  concerned  with  some  improvements  with 
respect  to  Ref.  2.  In  particular,  in  the  present  version  the 


aerodynamic  code  takes  into  account  the  wing  thickness  and 
the  real  fuselage  shape.  From  the  flight  mechanics  point  of 
view,  the  pitching  moment  at  zero  lift  has  now  been 
considered. 

The  whole  methodology  is  strongly  dependent  of  the 
optimisation  algorithm  both  from  the  computational 
performance  point  of  view  and  the  accuracy  of  the  obtained 
results.  Therefore,  a critical  analysis  of  the  available  methods 
has  been  carried  out:  this  is  described  in  the  next  section. 


2.  A CRITICAL  REVIEW  OF  OPTIMISATION 
METHODS 

The  main  limitation  to  aerodynamic  shape  optimisation  in 
a difficult  context  is  certainly  related  to  his  prohibitive 
computational  cost  due  to  complex  and/or  three  dimensional 
geometry,  and  to  multi-objective  problems  (multipoints 
and/or  multidisciplinary  optimisation).  Among  the  different 
methodologies  that  are  being  pursued  for  direct  optimisation, 
gradient-based  methods,  in  which  a specified  objective  is 
minimised,  are  often  employed.  These  procedures  require 
computations  of  both  objective  functional  and  its  derivatives 
respect  to  each  control  variable.  The  gradient  calculations  are 
usually  performed  using  black  box  finite  difference  methods 
which  are  simply  to  implement  but  require,  for  each  control 
variable,  at  least  a cost  evaluation.  For  practical  problems,  a 
large  number  of  control  variables  are  needed  to  take  care  of 
geometry  and  multidisciplinary  constraints.  Accordingly, 
when  the  functional  is  expensive  to  compute,  this  kind  of 
method  gives  poor  computational  performances. 

A first  possible  solution  is  to  use  relative  simple 
aerodynamic  models,  as  for  example,  potential  flows  (see  e.g. 
Refs.  4-6).  The  adjoint  approach,  based  on  control  theory  or 
variational  analysis,  is  an  interesting  alternative  procedure  for 
the  calculation  of  the  functional  derivatives,  which  makes  use 
of  a continuous  as  well  as  a discretised  differentiation  (see 
Newman  et  al.7  for  a review).  Actually,  it  allows  one  to 
greatly  limit  the  number  of  cost  evaluations,  and  thus,  can 
significantly  reduce  the  computational  cost.  For  the  above 
reasons,  in  the  last  years,  several  optimisation  computations 
for  complex  configurations  and  complex  flow  models  have 
been  performed  using  the  adjoint  method.  Usually,  this 
technique  is  combined  with  multigrid-like  algorithms  in  order 
to  reduce  the  dependency  between  cost  and  number  of  control 
variables  and/or  with  the  powerful  of  parallel  computing.  In 
this  way,  optimisations  of  realistic,  complete  aircraft 
geometries,  using  the  3D  Euler  equations  as  flow  model,  have 
been  performed811  as  well  as  2D  airfoil  shape  optimisations 
using  the  Reynolds  averaged  Navier-Stokes  equations 
coupled  with  turbulence  models  as  flow  model11'13. 

The  hand-code  exact  derivation  is  extremely  difficult  to 
perform  due  to  the  complexity  of  both  equations  and 
discretisation,  thus  automatic  differentiation  techniquesN  have 
been  used  for  the  evaluation  of  the  derivatives  in  a discrete 
adjoint  method13'15.  It  should  be  noted  that  methods  of  second- 
order,  .i.e.  methods  in  which  not  only  a gradient,  but  also,  a 
Hessian  must  be  evaluated,  can  be  found  in  literature  (see  e.g. 
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Ref.  16).  In  a context  of  multidisciplinary  design  problems, 
gradient  methods  minimise  a functional  constituted  by  a 
combination  of  different  criteria  with  different  weights. 
Alternative  methods,  which  do  not  require  the  use  of 
gradients,  can  be  also  considered,  in  particular  stochastic  ones 
as  the  genetic  algorithms  (GAs)  (see  Hajela17  for  a review). 
GAs,  which  have  found  significant  interest  for  applications  in 
the  aeronautical  field  (see  e.g.  Mosetti  and  Poloni18,  Crosslet 
and  Laananen19).  They  are  particularly  attractive  in  a complex 
context  because  of  their  robustness.  Indeed,  they  can 
overcome  the  limitations  of  gradient-based  methods,  i.e.  they 
can  operate  with  very  poor  regularity,  especially  with 
irregular  non- differentiable  functions  and  disjointed  feasible 
domains.  They  can  deal  with  a mix  of  different  objects  as 
continuous,  discrete,  integer  or  boolean  variables,  and,  they 
are  available  to  approach  a global  optimum  whereas  the 
gradient-based  methods  can  yield  only  local  optima. 
However,  the  major  drawback  of  GAs  is  that  they  require  a 
high  number  of  evaluations  of  the  objective  function,  and 
thus,  these  methods  are  poorly  efficient  in  a context  of 
expensive  cost  function.  In  particular,  once  the  nearness  of  the 
global  optimum  is  reached,  the  non  gradient  methods 
converge  very  poorly  whereas  this  condition  is  well-suited  for 
gradient-based  methods. 

To  partially  overcome  the  lack  of  computational 
efficiency,  parallel  computer  architectures  can  be  used  (see 
e.g.  Refs.  20-22)  drawing  profit  from  the  structure  of  GAs 
which  is  well  adapted  to  a fully  parallelisation  of  the 
algorithm.  Gradient-based  methods  and  GAs  have 
complementary  properties,  thus,  in  order  to  combine  their 
favourable  features,  an  interesting  approach  for  complex 
configurations  as  multidisciplinary  aeronautics  optimisation  is 
to  couple  the  two  methods.  Starting  from  a large  area  of 
configurations,  a GA  can  be  used  to  get  close  to  a particular 
configuration  which  is,  then,  optimised  by  a gradient-based 
method  to  reach  the  global  optimum.  This  kind  of  hybrid 
strategies  begins  to  be  studied  as  it  was  reported  in  very 
recent  works.23  24 

In  our  opinion  this  latter  approach  appears  to  be  the  most 
effective  for  the  problem  at  hand,  but  in  the  first  phase  of  our 
research  only  a gradient  method  has  been  applied.  Actually, 
this  method  has  been  used  to  obtain  the  simulation  results  for 
the  present  paper. 

3.  THE  MATHEMATICAL  MODEL 
3.1  Evaluation  of  aerodynamic  characteristics 

The  pressure  distribution  acting  on  the  airplane  was 
evaluated  by  means  of  a non  linear  panel  method,  based  on 
the  formulation  due  to  Morino;  the  features  of  this  method  are 
discussed  in  detail  in  Ref.  25. 

The  code  allows  for  wake  relaxation,  thus  yielding  non- 
loaded  stream  surfaces  as  required  by  the  correct  boundary 
condition.  Once  the  pressure  distribution  is  known,  the  lift 
distribution  over  wing,  fuselage  and  stabiliser  surfaces,  as  well 
as  the  induced  drag,  can  be  calculated. 


In  consideration  of  the  aerodynamic  code  chosen,  the 
analysis  is  naturally  confined  to  flows  at  low  angles  of  attack 
and  subsonic  Mach  numbers;  nevertheless,  it  can  provide 
useful  information  without  the  need  for  complex  and 
expensive  experimental  and/or  numerical  investigations. 

Aircraft  drag,  at  each  step  of  the  optimisation  procedure, 
was  estimated  by  using  the  component  building  method.  The 
induced  drag  of  the  configuration  was  predicted  by  means  of 
the  previously  described  code.  The  profile  drag  of  the  lifting 
surfaces  was  evaluated  by  means  of  standard  preliminary 
methodologies,  as  described  in  Refs.  26  and  27,  while  the 
drag  of  the  non  lifting  bodies  (fuselage  and  vertical  tail)  was 
evaluated  following  similar  methodologies.26,28  The  drag 
increment  caused  by  the  propulsive  system  installation  is  very 
sensitive  to  the  choice  of  a specific  configuration;  since  in  this 
phase  no  consideration  at  all  is  made  on  propulsion  system 
choice  and  position,  it  was  preferred  to  avoid  the  estimation  of 
its  effects.  Consequently,  the  total  drag  is  slightly 
underestimated,  even  if  this  simplification  does  not 
significantly  affect  the  optimisation  procedure.  Obviously, 
there  are  no  obstacles  in  taking  into  account  this  quantity 
where  more  detailed  studies  are  performed. 

3.2  The  numerical  optimisation  routine 

In  the  analysis  through  direct  numerical  optimisation,29  an 
aerodynamic  code  is  coupled  in  a loop  with  an  optimisation 
routine,  so  as  to  automatically  manage  the  values  of  the  design 
variables  - typically  concerning  geometry  modifications  - 
with  the  aim  of  minimising  a given  scalar  quantity  (objective 
function).  This  approach  is  extremely  flexible,  and  capable  of 
meeting  multidisciplinary  requirements.  The  choice  of  the 
design  variables  is  of  fundamental  importance,  as  they  define 
the  set  of  solutions  within  which  the  optimal  one  is  sought. 
Finally,  the  design  problem  formulation  requires  the  definition 
of  the  constraint  functions,  which  can  be  of  aerodynamic, 
flight  mechanics  or  geometrical  nature,  or  more  generally 
involve  any  quantity  that  may  be  computed  with  sufficient 
reliability. 

It  should  be  emphasised  that  a numerical  procedure  of  this 
kind  can  never  be  regarded  as  a completely  automated  tool; 
user’s  experience  and  external  control  on  the  process  are 
always  needed  to  obtain  the  best  results.  In  fact,  the  topology 
of  the  feasible  region  is  unknown  a priori. 

The  numerical  optimisation  routine  adopted  for  the  present 
study  is  the  CONMIN  code,30  used  extensively  in  the  fields  of 
aerodynamic  and  structural  optimisation.  CONMIN  uses  a 
gradient  method  for  the  search  algorithm  and  includes  three 
alternative  first  order  methods  for  the  calculation  of  the  vector 
search  direction.  When  all  the  constraints  are  satisfied  for  the 
current  values  of  the  design  variables,  the  steepest  descent 
method  is  used  for  the  first  iteration,  and  the  conjugate 
direction  method  for  subsequent  iterations.  When  some  of  the 
constraints  are  active  or  violated  the  feasible  directions 
method  is  used. 
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3.3  The  optimisation  procedure 

In  this  paper  the  main  objective  is  to  evaluate  the 
capabilities  of  the  proposed  methodology  to  obtain 
configurations  with  improved  performances  with  respect  to  a 
reference  one.  The  analysis  is  carried  out  upon  assuming  that 
the  following  design  parameters  have  been  arranged  for 
cruise:  payload,  velocity,  range,  cruise  height,  engine  specific 
fuel  consumption.  We  stress  that  the  different  configurations 
were  designed  only  for  achieving  a high-speed,  long  range 
cruise  mission  by  matching  the  airplane  design  to  cruise 
performance  design.  No  dynamic  analysis  nor  high  lift 
conditions  were  directly  considered.  However,  the  high  lift 
condition  is  implicitly  inserted  in  the  design  trough  a 
constraint  on  the  wing  loading,  as  discussed  in  the  following. 

It  has  been  supposed  that  payload  establishes  the  length 
and  diameter  of  the  fuselage,  which  have  been  considered 
known.  Thickness  ratio,  sweep  and  dihedral  angles  for  both 
wing  and  stabiliser,  and  the  ratio  between  the  wing  and  tail 
surfaces  have  been  assumed  as  design  parameters  and  have 
been  kept  fixed.  As  to  tail  volume  (imposed  here  by  means  of 
the  ratio  between  wing  and  tail  surfaces),  this  is  obviously 
imposed  by  mission  segments  different  from  cruise  (typically, 
take-off  and  landing).  The  sweep  angles  are  imposed  by  the 
Mach  number,  while  dihedral  angles  comes  from  lateral 
stability  considerations.  Finally,  as  to  thickness  ratio,  this  is  a 
function  of  both  cruise  and  high-lift  conditions.  Cruise 
imposes  a maximum  t/c  owing  to  drag  problems;  as  to  high  lift 
conditions,  it  can  be  easily  shown  that  the  function  cLmax 
versus  t/c  has  a maximum  which  depends  on  the  wing  section. 
This  entails  another  constraint  on  the  maximum  value  of  t/c . 
The  choice  of  some  of  the  above  design  parameters  may  have 
significant  effects  on  the  results;  accordingly,  the  influence  of 
a variation  in  thickness  ratio  and  tail  volume  coefficient  could 
be  investigated,  as  shown  in  Ref.  2. 

For  each  lifting  surface,  the  design  variables  used  in  the 
optimisation  process  are  aspect  ratio,  taper  ratio  and  twist. 
Moreover,  wing  loading  and  stabiliser  position  (horizontal 
stagger  and  vertical  gap)  are  also  design  variables.  However, 
wing  loading  is  a variable  which  is  typically  dependent  upon 
flight  phases  other  than  cruise  and,  therefore,  it  is  not  possible 
to  optimise  this  variable  without  a constraint  which  forbids  it 
to  assume  unrealistic  values.  Moreover,  stabiliser  position  is 
constrained  by  fuselage  geometry,  since  its  length  and 
diameter  have  been  fixed. 

It  is  well  known  that  one  of  the  crucial  aspects  of  the 
optimisation  procedure  is  constituted  by  the  choice  of  the 
object  function.  In  the  present  analysis  various  considerations 
lead  us  to  the  conclusion  that  the  take-off  weight  is  a 
reasonable  choice,  even  if  it  should  be  clear  that  the  essence 
of  the  proposed  methodology  would  not  be  affected  by  a 
different  objective  function.  The  main  reason  is  that  the 
aircraft  take-off  weight  appears  to  offer  a quite  reasonable 
compromise  between  the  need  to  obtain  solutions  that  are 
sufficiently  accurate  to  give  significant  results  and  the  need  to 
prevent  such  solutions  from  being  excessively  complicated,  as 
it  would  happen  by  using,  for  example,  Direct  Operating 


Costs  as  an  objective  function.31 

The  objective  function  is  considered  to  be  the  sum  of  four 
terms:  wing  weight,  stabiliser  weight,  fuel  weight  and  residual 
weight.  The  latter  comprises  fuselage,  payload,  vertical  tail, 
engines,  gears  and  various  systems;  their  weights  and  centre 
of  gravity  positions  are  assumed  to  be  known,  and  have  been 
established  by  statistical  data. 

The  optimisation  process  includes  the  following  steps: 

1.  On  the  basis  of  the  above  design  parameters  and  with  a 
fixed  set  of  design  variables,  a first  estimation  of  WD  and 
W/ is  given  and,  as  a function  of  W0  and  Wf,  an  estimation 
of  Ww  and  Ws  is  obtained  by  using  formulae  derived  from 
statistical  data.32  Accordingly,  wing  and  stabiliser 
planform  is  fixed  as  a function  of  the  distance  between 
wing  and  stabiliser  aerodynamic  centres.  Centre  of  gravity 
of  wing,  stabiliser  and  fuel  (supposed  to  be  concentrated 
in  the  wing)  are  also  evaluated.  In  the  following  examples 
the  positions  of  the  centres  of  gravity  of  various  elements 
are  obtained  by  means  of  simplified  formulae.33 

2.  The  angle  of  attack  and  elevator  angle  are  computed  to 
trim  the  aircraft  in  straight,  horizontal  flight,  by  using  the 
following  equations 

2 W , 

CL,r,m  = „ c v2  = C*-Lo,S=0  + CLaCl  + CLS5  (l) 

Pc^w*c 

C M ~ 0 — CM  |a_Q  ^_Q  +cMa  & J<rCM  (2) 

The  aerodynamic  coefficients  in  Eqs.  (1)  and  (2)  are 
computed  with  an  iterative  procedure,  by  using  the 
previously  described  non-linear  aerodynamic  code.  The 
value  of  cDtrim  is  then  evaluated  with  the  previously 
described  aerodynamics  procedure  and  vehicle  efficiency 
is  obtained  by  E=cLtrJcDtrim. 

3.  By  means  of  the  computed  values  of  cLa  and  cMa  the  static 
longitudinal  stability  of  the  configuration  is  checked.  The 
code  constrains  the  difference  between  centre  of  gravity 
and  aerodynamic  centre  positions  (non-dimentionalised  by 
cma ) to  be  included  between  two  imposed  bounds.  In  the 
examples  the  inequality  -0.05  < (h-hn)  < -0.2  was  used. 
When  the  above  constraint  is  violated,  the  code 
repetitively  comes  into  play  to  modifying  the  value  of  lt , 
by  moving  the  wing  with  respect  to  fuselage,  until  the 
bounds  are  met. 

4.  At  this  point,  the  load  distribution  on  lifting  surfaces  is 
known,  therefore  it  is  possible  to  give  a new  estimation  of 
wing  weight  by  using  the  methodology  of  Torenbeek3 
which  binds  wing  weight  to  lift  distribution  along  the 
wing  span  and  mechanics  characteristics  of  the  employed 
materials  (see  Appendix).  By  this  means  it  is  possible  to 
accurately  estimate  wing  weight. 

Fuel  weight  is  obtained  from  Breguet’s  formula: 
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Wf=W0 


1 - exp 


f RSfcS 

VCE 


(3) 


and  the  quantity  W ~Ww+Ws+Wj+Wres  is  evaluated.  If 
W ±W0  (except  for  numerically  insignificant 
differences),  a new  estimation  of  WQ  is  given  by  using,  for 
example,  the  updating  formula  WQ  = (w  + W0)/2 . 


5.  The  points  1-4  are  repeated  until  the  condition  W = W0  is 
met.  It  should  be  noted,  however,  that  Ww  at  point  1 is 
now  evaluated  by  Torenbeek’s  approach. 

6.  The  optimisation  routine  updates  the  design  variables  and 
the  procedure  described  at  points  1-5  is  repeated  until  the 
minimum  weight  configuration  is  obtained. 

We  stress  the  fact  that,  since  the  wing  weight  is  evaluated 
by  means  an  approach  based  on  the  real  span  wise  lift 
distribution,  (and  not  simply  on  statistical  data)  it  is  possible 
to  analyse  canard  configuration  as  well.  Nevertheless,  in  this 
case  the  estimations  based  upon  a statistical  approach  are 
defective  due  to  the  small  number  of  available  data. 


given  by  the  constraint.  On  the  other  hand,  the  optimisation  of 
tail  variables  turned  out  to  have  a negligible  effect,  and  the  tail 
position  reached,  in  all  the  cases  presented,  the  upper  values 
of  the  constraints.  The  most  important  results  of  the 
optimisation  process  have  been  summarised  in  Table  2a, 
where  aircraft  geometry  is  defined,  and  Table  2b,  where  trim 
conditions  are  shown. 


Wing 

Tail 

s 

23.46 

5.87 

b 

15.14 

5.97 

cr 

1.80 

1.31 

AR 

9.76 

6.07 

A 

0.72 

0.5 

e 

-2.64° 

0° 

h 

5.46 

10.9 

T 

2.00 

i 

0( 

3 

Ww 

5461 

Ws 

514 

Wf 

5686 

W0 

44652 

Table  2a  - Aircraft  geometry  and  weights 


4.  THE  OPTIMISATION  PROCEDURE:  A CASE 
STUDY 

In  this  section  the  outlined  methodology  has  been  applied, 
as  an  example,  to  a typical  light  transport  aircraft.  Table  1 
shows  the  values  of  the  fixed  parameters  used  in  the 
optimisation  procedure. 


Wp  xce 

M 

R 

Sft 

H 

33000  40 

0.4 

2100 

0.48 

6100 

Lfus  fyus 

A, 

A 

'h 

rf 

11.2  1.7 

0° 

4° 

4° 

9° 

Table  1 - Fixed  parameters  used  in  the  optimisation 
procedure 


In  addition  to  the  data  of  Table  1,  the  following  design 
parameters  have  been  fixed:  (t/c)  from  0.18  (at  the  root)  to 
0.12  (at  the  tip),  both  for  wing  and  tail.  Moreover  a ratio 
S/SH= 0.25  has  been  fixed. 

The  design  variables  are  aspect  ratio,  taper  ratio  and  twist 
of  wing  and  tail,  wing  loading  and  tail  position.  The 
constraints  are  as  follows:  1)  maximum  value  of  wing  loading 
WS=1903  N/m2;  2)  the  tail  vertical  position  could  only  vary 
between  the  maximum  height  of  the  fuselage  and  two  meters 
above  that  position;  3)  tail  horizontal  position  constrained  by 
fuselage  length.  Note  that  a fully  rotating  tail  has  been 
assumed.  It  is  interesting  to  note  that  the  objective  function 
was  highly  sensitive  to  wing  loading  which  reached  the  value 


a 

1.52° 

S 

1.37° 

Cl 

0.361 

Cd 

0.0218 

E 

16.6 

Table  2b  - Trim  conditions 


The  geometry  of  the  optimised  configuration  is  shown  in 
Fig.  1 along  with  the  mesh  used  in  the  computational  analysis 
(for  the  body  representation,  about  2700  panels  have  been 
used).  Typical  computational  times  required  to  obtain  an 
optimised  configuration  are  of  the  order  of  5h  on  a Pentium 
III  Xion,  500  MHz  and  512  MB  RAM. 

As  far  as  drag  analysis  is  concerned,  wing  loading  being 
practically  fixed  by  the  constraint,  it  can  be  seen  that  wetted 
area  is  practically  a constant  (it  is  clear  that,  for  a correct 
optimisation  process,  the  final  take-off  weight  cannot  be 
significantly  different  from  the  initial  estimation),  hence 
friction  drag  does  not  vary  to  any  remarkable  extent  and  the 
optimisation  process  tends  almost  exclusively  to  decrease  the 
induced  drag.  Accordingly,  the  optimised  configuration  is  the 
best  compromise  between  two  opposite  standpoints:  an 
increase  in  wing  span  to  decrease  induced  drag  (thereby 
decreasing  Wf  and  W0)  or,  conversely,  a reduction  in  wing 
span  to  reduce  structural  stress  and,  as  a result,  to  diminish  Ww 
and  W0. 


Fig.  1 - Reference  optimised  configuration 


5.  A COMPARISON  OF  DIFFERENT  CHOICE  OF 
THE  GEOMETRICAL  PARAMETERS 

Although  the  above  parameters  are  constrained,  from  a 
practical  point  of  view,  by  mission  segments  different  from 
cruise  (especially  take  off  and  landing),  it  is  interesting  to 
obtain  quick  estimates  of  the  order  of  magnitude  of  the 
improvements  which  would  be  related  to  a relaxation  of  these 
constraints.  This  in  turn  gives  important  information  about 
the  trade-off  between  an  increase  in  costs  related  to  more 
sophisticated  solutions  and  an  increase  in  the  overall 
configuration  performance. 

The  optimised  configuration  outlined  in  the  previous 
section  results  from  a number  of  different  geometrical 
constraints.  It  is  interesting  to  investigate  the  effect  of  the 
parameters  which  are  not  directly  controlled  by  the  optimiser. 
In  particular,  we  analyse  the  differences  in  the  configuration 
due  to  a different  choice  of  the  wing  thickness  ratio  and  to  a 
different  wing  loading. 

In  Table  3 the  optimised  configuration  obtained  with  a 
greater  wing  loading  constraint  is  compared  to  the  reference 
one. 

A significant  reduction  in  wing  weight  with  an  increase  in 
wing  loading  is  apparent,  along  with  a less  important  effect 
on  fuel  consumption.  Note  that  the  increase  in  wing  loading 
causes  a decrease  in  wing  span  and  an  slight  increase  in 
aerodynamic  efficiency,  thus  leading  to  both  a reduction  in 
wing  weight  and  lower  fuel  consumption. 


Value 

WS=  1903  (Ref) 

W5=2093 

S„ 

23.46 

20.93 

K 

15.14 

14.04 

AK 

9.76 

9.41 

K 

0.72 

0.72 

eK 

-2.64 

-4.64 

bs 

5.97 

5.64 

K 

0.5 

0.5 

es 

0 

0 

i, 

5.46 

5.14 

cL 

0.361 

0.396 

Cd 

0.0218 

0.0237 

E 

16.6 

16.7 

WK 

5461 

4752 

W5 

514 

461 

Wf 

5686 

5549 

W0 

44652 

43813 

(Wo-  woref)/  woref 

- 

-1.88% 

Table  3 - Effect  of  a different  wing  loading 


In  Table  4 the  optimised  configuration  obtained  with  a 
reduced  value  of  thickness  ratio  at  wing  root  is  compared  to 
the  reference  one.  In  particular,  the  reference  value, 
corresponding  to  t/c- 0.18  at  wing  root,  is  modified  in 
f/c=0.16. 
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Value 

t/cr00i — 0.18 

(Ref.) 

t/c  root—  0.16 

Sw 

23.46 

23.54 

bK 

15.14 

15.27 

ark 

9.76 

9.90 

K 

0.72 

0.71 

eK 

-2.64 

-2.68 

bs 

5.97 

5.98 

A 

0.5 

0.5 

es 

0 

0 

h 

5.46 

5.34 

cl 

0.361 

0.360 

CD 

0.0218 

0.0220 

E 

16.6 

16.4 

w„ 

5461 

5505 

514 

524 

wf 

5686 

5763 

44652 

44806 

(W„-  Wnref)l  wnref 

- 

+0.34% 

Table  4 - Effect  of  a different  thickness  ratio 


It  is  clear  that  weight  increases  as  t/c  at  the  root  decreases. 
Upon  decreasing  t/c  at  the  root,  performance  decreases  due  to 
an  increase  in  wing  weight.  It  should  be  noted,  however,  that 
the  higher  structural  efficiency  is  not  “used”  to  increase  wing 
span,  which,  indeed,  shows  a slight  decrease.  When  t/c  at  the 
root  is  increased  from  0.16  up  to  0.18,  the  higher  structural 
efficiency  is  “used”  to  increase  wing  span,  thus  reducing  the 
induced  drag,  with  a consequently  lower  fuel  consumption, 
but  without  important  differences  in  wing  weight. 

6.  A COMPARISON  OF  DIFFERENT  CHOICE  OF 
THE  TECHNOLOGICAL  PARAMETERS 

In  this  section,  the  outlined  methodology  has  been  applied 
to  stress  the  importance  of  some  of  the  most  significant 
technological  parameters  over  the  final  aircraft  configuration. 

6.1  Reduction  in  friction  drag 

It  is  known  that  the  reduction  of  friction  drag  plays  a 
relevant  role  in  the  current  aerodynamic  research.  In  what 
follows,  we  suppose  that  the  friction  drag  may  be  reduced  by 
10%  with  respect  to  the  nominal  case.  As  a first  example,  this 
reduction  has  been  applied  to  the  optimised  configuration  of 
Table  2,  which  results  in  a fuel  reduction  of  370  N,  which  in 
turn  corresponds  to  a reduction  of  0.8%  with  respect  to  the 
take  off  weight. 

It  is  interesting  to  investigate  if  the  optimisation  process 
would  be  different  in  this  case:  accordingly  the  whole 
optimisation  procedure  has  been  carried  out  with  the  new 
value  of  the  friction  coefficient  and  the  results  are  shown  in 
Table  5. 

The  optimised  configuration  obtained  considering  a 
smaller  friction  drag  shows  a reduction  in  the  take  off  weight 
of  1.33%,  significantly  higher  (+0.53%)  than  the 


corresponding  result  obtained  considering  the  reduction  in 
friction  drag  applied  to  the  reference  configuration. 


Value 

Reference 

cf-0.9cfivr 

SK 

23.46 

23.15 

bw 

15.14 

14.64 

AR  '»• 

9.76 

9.26 

A, 

0.72 

0.74 

0, 

-2.64 

-4.37 

bs 

5.97 

5.93 

A 

0.5 

0.5 

es 

0 

0 

i, 

5.46 

5.38 

cl 

0.361 

0.360 

Cd 

0.0218 

0.0206 

E 

16.6 

17.5 

ww 

5461 

5166 

514 

505 

wf 

5686 

5330 

W0 

44652 

44057 

(W„-  W„rrf)/  wnrtf 

- 

-1.33% 

Table  5 - Effect  of  a different  friction  drag 


The  new  configuration  shows  an  increased  aerodynamic 
efficiency,  even  with  a reduced  wing  span;  therefore,  both  the 
wing  and  fuel  weights  are  reduced.  The  lift  coefficient 
remains  practically  constant,  therefore  the  aerodynamic  flow 
is  practically  the  same,  while  the  geometric  configuration 
appears  slightly  different.  Indeed,  the  new  configuration  has  a 
reduced  wing  surface  with  a reduction  both  in  wing  span 
(-3.3%)  and  aspect  ratio  (-5.1%),  a slightly  higher  taper  ratio 
and  an  increased  twist. 

6.2  Reduction  in  specific  fuel  consumption 

Another  fundamental  research  area  in  the  aeronautical 
field  relates  to  the  increase  in  engine  efficiency.  Accordingly, 
in  what  follows  we  suppose  that  the  specific  fuel  consumption 
may  be  reduced  by  10%  with  respect  to  the  nominal  case.  As 
a first  example  this  reduction  has  been  applied  to  the 
optimised  configuration  of  Table  2,  which  results  in  a fuel 
reduction  of  569  N,  which  corresponds  to  a reduction  of  1.3% 
with  respect  to  the  take  off  weight.  Obviously,  in  this  case  the 
reduction  in  fuel  weight  is  very  close  to  the  corresponding 
reduction  in  specific  fuel  consumption. 

As  in  the  previous  example,  the  whole  optimisation 
procedure  has  been  carried  out  with  the  new  value  of  the  new 
value  of  the  specific  fuel  consumption  and  the  results  are 
shown  in  Table  6. 

The  optimised  configuration  obtained  considering  a 
smaller  specific  fuel  consumption  shows  a reduction  in  the 
take  off  weight  of  1.98%,  significantly  higher  (+0.68%)  than 
the  corresponding  result  obtained  considering  the  reduction  in 
specific  fuel  consumption  applied  to  the  reference 
configuration.  The  reduction  in  specific  fuel  consumption 
directly  reduces  the  fuel  weight,  with  a consequent  reduction 
in  the  wing  weight,  mainly  caused  by  the  reduction  in 


required  total  lift. 


Value 

Reference 

S',=0.9Str„ 

5, 

23.46 

23.00 

ft 

15.14 

14.59 

ark 

9.76 

9.26 

ft 

0.72 

0.74 

ft 

-2.64 

-4.55 

ft 

5.97 

5.91 

ft 

0.5 

0.5 

ft 

0 

0 

h 

5.46 

5.38 

cl 

0.361 

0.360 

Cd 

0.0218 

0.0220 

E 

16.6 

16.4 

WK 

5461 

5135 

w. 

514 

501 

Wf 

5686 

5119 

W0 

44652 

43767 

(W0-Woref)IWoref 

-1.98% 

Table  6 - Effect  of  a different  SfC 


The  new  configuration  shows  an  aerodynamic  behaviour 
practically  unchanged  with  respect  to  the  reference 
configuration:  the  lift  coefficient,  the  drag  coefficient  and  the 
efficiency  remain  practically  constant,  therefore  the 
aerodynamic  flow  is  almost  the  same,  while  the  geometric 
configuration  appears  to  be  slightly  different.  Indeed,  the  new 
configuration  has  a reduced  wing  surface  with  a reduction 
both  in  wing  span  (-3.3%)  and  aspect  ratio  (-5.1%),  a slightly 
higher  taper  ratio  and  an  increased  twist. 

6.3  Increase  in  structural  efficiency 

As  a final  example,  we  investigate  the  effect  of  an  increase 
of  10%  of  the  admissible  normal  stress  level  of  the  material 
over  the  optimised  configuration.  Upon  applying  this 
modified  value  to  the  optimised  configuration  of  Table  2,  a 
reduction  of  wing  weight  of  277  N is  obtained  along  with  a 
reduction  of  35  N of  fuel  weight,  which  corresponds  to  a 
reduction  of  0.6%  with  respect  to  the  take  off  weight. 

The  whole  optimisation  procedure  has  been  carried  out 
with  the  new  value  of  oa  and  the  results  are  shown  in  Table  7. 

The  optimised  configuration  obtained  considering  a higher 
structural  efficiency  shows  a reduction  in  the  take  off  weight 
of  0.93%,  higher  (+0.33%)  than  the  corresponding  result 
obtained  considering  the  increase  in  structural  efficiency 
applied  to  the  reference  configuration,  but  not  so  important  as 
in  previous  analysed  cases. 

The  new  configuration  shows  an  aerodynamic  behaviour 
practically  unchanged  with  respect  to  the  reference 
configuration:  the  lift  coefficient,  the  drag  coefficient  and  the 
efficiency  remain  practically  constant,  therefore  the 
aerodynamic  flow  is  almost  the  same.  The  geometric 
configuration  appears  substantially  unchanged,  with  an 
increased  twist  as  the  only  significant  modification. 


Value 

Reference 

Off  —1-1  0({  ri,f 

5, 

23.46 

23.24 

K 

15.14 

15.08 

ARW 

9.76 

9.79 

ft.- 

0.72 

0.72 

ft. 

-2.64 

-4.49 

ft 

5.97 

5.94 

ft 

0.5 

0.5 

0s 

0 

0 

l, 

5.46 

5.29 

cL 

0.361 

0.360 

Cd 

0.0218 

0.0218 

E 

16.6 

16.5 

WK 

5461 

5024 

514 

508 

Wf 

5686 

5656 

Wo 

44652 

44235 

(W„-  Woref)l  wmef 

- 

-0.93% 

Table  7 - Effect  of  a different  oa 


Therefore,  the  reduction  in  the  take  off  weight  is  mainly 
produced  by  the  direct  effect  of  the  increased  structural 
efficiency  on  the  wing  weight  and  a consequent  limited  effect 
on  the  fuel  weight. 

7.  CONCLUDING  REMARKS 

A methodology  for  the  analysis  and  comparison  of  the 
performance  of  different  configurations  in  the  phase  of 
aircraft  conceptual  design  has  been  developed.  In  the 
analysis,  a scalar  objective  function,  the  take-off  weight,  is 
minimised  by  means  of  a numerical  optimisation  technique, 
which  takes  into  account  the  high  number  of  geometrical 
parameters  and  the  flight  mechanics  requirements  involved  in 
the  problem.  The  approach  proved  to  be  extremely  flexible, 
and  capable  of  really  meeting  multidisciplinary  requirements. 

The  obtained  results  highlight  the  order  of  magnitude  of 
the  improvements  of  the  objective  function  due  to  an  increase 
in  technological  characteristics.  Moreover,  it  is  clear  that 
significantly  higher  improvements  in  the  performances  can  be 
obtained  if  the  configuration  is  developed  by  taking  into 
account  the  above  technological  improvements  from  the 
initial  phase  of  the  project.  This  means  that  the  optimised 
final  configuration  is  different  from  the  reference  one  (which 
in  turn  has  been  optimised  with  different  values  of  these 
parameters). 

We  also  note  that  the  obtained  results  are  “close”  to  the 
reference  configuration.  This  is  due  to  the  fact  that  a gradient 
based  optimiser  may  only  give  local  minima.  Actually,  better 
results  may  be  obtained  by  combining  a global  (GA  based) 
optimiser,  with  a local  gradient  based  approach.  This  matter  is 
currently  under  development. 

It  is  clear  that  the  obtained  results  are  not  definitive,  since 
a significant  number  of  different  configurations  must  be 
considered  and  several  points  need  to  be  verified  with  a more 
accurate  analysis.  Among  them  the  most  important  one  is  the 
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high  lift  condition  which  was  not  considered  directly  in  the 
optimisation  process,  even  if  we  assumed  that  it  imposes  a 
wing  loading  constraint.  Obviously,  this  constraint  must  be 
accurately  verified  with  an  a posteriori  analysis. 


8.  APPENDIX 

According  to  Torenbeek’s  approach,3  basic  wing  weight  is 
considered  as  the  sum  of  three  contributions:  1)  the  weight  of 
material  required  to  resist  bending,  WBL  2)  the  weight  of 
material  required  to  resist  shear  forces,  WSl  and  3)  the  weight 
of  wing  ribs,  W^.  All  three  contributions  are  calculated  as 
follows. 


8.3  Weight  of  wing  ribs 

Since  a rational  derivation  for  the  weight  of  ribs  is  not 
feasible,  the  following  (statistical)  formulation  is  used: 

Wfib  :=  ^rPM  8^  * [tref  ^ r ^ h V 2] 

where  kr  =0.5xl(T3  and  trej  =1.0  m. 

Finally,  basic  wing  weight  estimation  is  the  sum  of  the 
above  three  contributions,  hence: 

Ww=WBL+WSL+Wrib 


REFERENCES 


8.1  Weight  due  to  bending  moment 

The  bending  moment  due  to  lift  at  a generic  (y)  wing 
station  is  given  by: 

MBL(y)=  [ — y qci(y')c(y')dy' 

Jy  COSAea 

and  c/(y')  is  evaluated  by  the  aerodynamic  code. 

The  basic  weight  of  two  skin  panels  (upper  and  lower)  can 
be  estimated  by  the  formula: 


with  A(y)  given  by: 


A(y)  = 


MBL(y) 

n,‘(y)Ga 


where  Ga  is  the  given  admissible  normal  stress  level,  t(y)  is 
the  maximum  thickness  of  the  profile  section  and  rjt  is  an 
efficiency  factor  which  may  be  obtained  from  a drawing  of 
the  wing  cross  section  or  by  simplified  expressions.3 

8.2  Weight  due  to  shear  forces 

The  weight  due  to  shear  forces  can  be  expressed  as 

WSL=^£fi-MBL(y=  0) 


where  T is  the  mean  value  of  shear  stress  which  has  been 
assumed  constant,  f can  be  estimated  by  means  of  the  simple 
relationship: 


t = 0.5cra 


1.  Lombardi,  G.  and  Vicini,  A.,  “Induced  Drag  Prediction 
for  Wing-Tail  and  Canard  Configurations  through 
Numerical  Optimization,’’  The  Aeronautical  Journal, 
Vol.  98,  No.  976  June-July  1994,  pp  199-206. 

2.  Lombardi  G.,  Mengali  G.  “A  Methodology  for  the 
Preliminary  Analysis  and  Comparison  of  Wing-Tail 
and  Canard  Configurations”,  The  Aeronautical  Journal, 
Vol.  101,  No.  1004,  April  1997,  pp.169-177. 

3.  Torenbeek,  “Development  and  Application  of  a 
Comprehensive,  Design-Sensitive  Weight  Prediction 
Method  for  Wing  Structures  of  Transport  Category 
Aircraft”,  TU  Delft,  Report  LR-693,  September  1992. 

4.  Selmin  V.,  “Optimization  of  Transonic  Arfoils”,  Notes 
on  Numerical  Fluid  Mechanics  Vol.  61,  Optimum 
Aerodynamic  Design  and  Parallel  Navier-Stokes 
Computations,  EC  ARP  - European  Computational 
Aerodynamics  Researsh  Project,  Periaux  et  al  Editors, 
1998,  pp.  60-70. 

5.  Fomassier  L.,  “Optimization  of  the  Wing-Pylon- 
Nacelle  Testcase  TE5  by  HISSS-D,  a Panel  Method- 
Based  Design  Tool”,  Notes  on  Numerical  Fluid 
Mechanics  Vol.  61,  Optimum  Aerodynamic  Design 
and  Parallel  Navier-Stokes  Computations,  ECARP  - 
European  Computational  Aerodynamics  Researsh 
Project,  Periaux  et  al  Editors,  1998,  pp.  99-108. 

6.  Labrujere  Th.E.,  “Residual-correction  type  and  related 
computational  methods  for  aerodynamic  design”, 
AGARD  Report  R-803  Optimum  Design  Methods  for 
Aerodynamics,  1994,  pp.  4-1/5-31. 

7.  Newman  III  J.C.,  Taylor  III  A.C.,  Barnwell  R.W., 
Newman  P.A.,  Hou  G.J.-W.,  “Overview  of  Sensitivity 
Analysis  and  Shape  Optimization  for  Complex 
Aerodynamic  Configurations”,  Journal  of  Aircraft,  Vol. 
36,  No.  1,  1999,  pp.  87-96. 

8.  Reuther  J.J.,  Jameson  A.,  Alonso  J.J.,  Rimlinger  M.J., 
Saunders  D.,  “Constrained  Multipoint  Aerodynamic 
Shape  Optimization  Using  an  Adjoint  Formulation  and 
Parallel  Computers”,  Journal  of  Aircraft,  Vol.  36,  No. 
1,  1999,  pp.  51-74. 

9.  Elliott  J.,  Peraire  J.,  “Practical  Three-Dimensional 
Aerodynamic  Design  and  Optimization  Using 


6-10 


Unstructured  Meshes”,  AIAA  Journal,  Vol.  35,  No.  9, 
1997,  pp.  1479-1485. 

10.  Marco  N.,  Dervieux  A.,  “’Multilevel  parametrization 
for  aerodynamical  optimization  of  3D  shapes”,  Finite 
Elements  in  Analysis  and  Design,  Vol.  26,  1997,  pp. 
259-277. 

11.  Newman  III  J.C.,  Taylor  III  A.C.,  Barnwell  R.W., 
“Aerodynamic  Shape  Sensitivity  Analysis  and  Design 
Optimization  of  Complex  Configurations  Using 
Unstructured  Grids”,  AIAA  Paper  97-2275,  1997. 

12.  Anderson  W.K.,  Bonhaus  D.L.,  “Airfoil  Design  on 
Unstructured  Grids  for  Turbulent  Flows”,  AIAA 
Journal,  Vol.  37,  No.  2,1999,  pp.  185-191. 

13.  Mohammadi  B.,  “Optimal  Shape  Design,  Reverse 
Mode  of  Automatic  Differentiation  and  Turbulence”, 
AIAA  Paper  97-0099,  1997. 

14.  Bischof  C.H.,  Carle  A.,  Corliss  G.F.,  Griewank  A., 
Hovland  P.,  “ADIFOR:  Generating  derivative  code 
from  FORTRAN  programs”,  Scientific  Programming, 
Vol.  1,  1992,  pp.  11-29. 

15.  Male  J.-M.,  Rostaing-Schmidt  N.,  Marco  N., 

“Automatic  Differentiation:  an  Application  to 

Optimum  Shape  Design  in  Aeronautics”, 
ECCOMAS'96  minisymposium  on  Automatic 
Differentiation  of  Programs  and  Large  Scale 
Optimization,  1996,  pp  87-91. 

16.  Jou  W.H.,  Huffman  W.P.,  Young  D.P.,  Melvin  R.G., 
Bieterman  M.B.,  Hilmes  C.L.,  Johnson  F.T.,  “Practical 
considerations  in  aerodynamic  design  optimization” 
AIAA-95-1730-CP,  1995,  pp.  950-960. 

17.  Hajela  P.,  “Nongradient  Methods  in  Multidisciplinary 
Design  Optimization  - Status  and  Potential”,  Journal  of 
Aircraft,  Vol.  36,  No.  1,  1999,  pp.  255-265. 

18.  Mosetti  G.,  Poloni  C.,  “Aerodymanic  shape 
optimization  by  means  of  a genetic  algorithm”,  Proc.  of 
the  5th  Symp.  on  Computational  Fluid  Dymanics, 
Sendai,  Vol.  II,  1993,  pp.  279-284. 

19.  Crosslet  W.A.,  Laananen  D.H.,  “Conceptual  Design  of 
Helicopters  via  Genetic  Algorithm”,  Journal  of 
Aircraft,  Vol.  33,  No.  6,  1996,  pp.  1062-1070. 

20.  Obayashi  S.,  Yamaguchi  Y.,  Nakamura  T., 

“Multi  objective  Genetic  Algorithm  for 

Multidisciplinary  Design  of  Transonic  Wing 
Planform”,  Journal  of  Aircraft,  Vol.  34,  No.  5,  1997, 
pp.  690-693. 

21.  Marco  N.,  Lanteri  S.,  Desideri  J.A.,  Periaux  J.,  “A 
Parallel  Genetic  Algorithm  for  Multi-Objective 
Optimization  in  Computational  Fluid  Dynamics”, 
Evolutionary  Algorithms  in  Engineering  and  Computer 
Science,  Miettinen  et.  al.  Eds,  J.  Wiley  & Sons,  1999, 
pp.  445-456. 

22.  Lee  J.,  Hajela  P.,  “Parallel  Genetic  Algorithm 
Implementation  in  Multidisciplinary  Rotor  Blade 
Design”,  Journal  of  Aircraft,  Vol.  33,  No.  5,  1996,  pp. 
962-969. 

23.  Dulikravich  G.S.,  Martin  T.J.,  Dennis  B.H.,  Foster 
N.F.,  “Multidisciplinary  Hybrid  Constrained  GA 
Optimization”,  Evolutionary  Algorithms  in 


Engineering  and  Computer  Science,  Miettinen  et.  al. 
Eds,  J.  Wiley  & Sons,  1999,  pp.  233-259. 

24.  Vicini  A.,  Quagliarella  D.,  “Airfoil  and  Wing  Design 
Through  Hybrid  Optimization  Strategies”,  AIAA 
Journal,  Vol.  37,  No.  5,  1999,  pp.  634-641. 

25.  Morino,  L.,  Kuo,  C.C.,  “Subsonic  Potential 
Aerodynamics  for  Complex  Configurations:  a General 
Theory”,  AIAA  Journal,  Vol.  12,  No.  2,  1974,  pp.  191- 
197. 

26.  Anonymous  “Usaf  Stability  and  Control  DATCOM”, 
item  4.1.5. L,  April  1978. 

27.  Anonymous  “Usaf  Stability  and  Control  DATCOM”, 
item  4.I.5.2.,  April  1978. 

28.  Anonymous  “Usaf  Stability  and  Control  DATCOM”, 
item  4.2.3. 1-A.,  April  1978. 

29.  Vanderplaats,  G.N.,  “Numerical  Optimization 

Techniques  for  Engineering  Design:  with 

Applications”,  Mac  Graw-Hill  series  in  Mechanical 
Engineering,  1984. 

30.  Vanderplaats,  G.N.,  “CONMIN,  a Fortran  program  for 
constrained  function  minimization”,  NASA  TM  X-62- 
282,  1973. 

31.  Knowles  and  Martinez-Val  “Conventional  and 
Unconventional  Configurations  for  Ultra-High 
Capacity  Aircraft”  ICAS  paper  94-1.3.2,  1994,  pp.656- 
664. 

32.  Raymer  D.P.,  “Aircraft  Design:  a Conceptual 

Approach”,  AIAA  Education  Series,  Washington  DC, 
1989. 

33.  Roskam,  J.  “Airplane  Design”,  Roskam  Aviation  and 
Engineering  Corporation,  Ottawa,  1990,  Part  VI  and 
VII. 


